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HIGHLIGHTS 


►  Three  start-up  procedures  for  an  anode-supported  planar  SOFC  are  investigated. 

►  The  largest  absolute  temperature-gradient  is  exhibited  in  the  early  stage. 

►  The  anode-recycling  start-up  procedure  significantly  reduces  the  start-up  time. 

►  The  fixed-temperature-difference  procedure  may  reduce  the  temperature-gradient. 

►  The  fixed-temperature-difference  procedure  also  reduces  the  start-up  time. 
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Three  start-up  procedures  for  an  anode-supported  planar  SOFC  are  proposed  and  investigated  numeri¬ 
cally  in  the  present  study.  The  first  is  to  introduce  the  inlet  fuel  at  the  operation  temperature  after  the 
heat-up  process  is  completed.  The  second  is  to  incorporate  the  anode-recycling  mechanism  into  the 
start-up  process.  The  third  is  to  fix  the  difference  between  the  inlet-fuel  temperature  and  the  cell 
minimum  temperature.  The  numerical  results  obtained  from  the  present  study  show  that  the  effective 
maximum  absolute  temperature-gradient  is  exhibited  in  the  early  stage  of  the  start-up  process.  For  the 
present  investigated  SOFC  configuration,  the  required  start-up  time  for  the  case  using  methane  is  3.2- 
fold  longer  than  that  using  hydrogen.  The  effective  maximum  absolute  temperature-gradient  for  the 
case  using  hydrogen  is  2.2-fold  larger  than  that  utilizing  methane.  The  endothermic  internal  reforming 
reaction  of  methane  has  a  positive  effect  on  the  accommodation  of  the  temperature  uniformity  during 
the  start-up  process.  The  anode-recycling  mechanism  significantly  reduces  the  start-up  time  For  the 
fixed-temperature-difference  start-up  procedure,  a  properly  selected  temperature  difference  may  lead  to 
a  smaller  effective  maximum  absolute  temperature-gradient  in  the  early  stage  and  a  shorter  start-up 
time  by  accelerating  the  start-up  pace  in  the  later  stage. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  solid  oxide  fuel  cell  (SOFC),  which  directly  converts  the 
chemical  energy  into  electrical  power  via  the  electrochemical 
reaction,  is  one  of  the  most  attractive  candidates  for  the  next 
generation  of  power  sources  [1].  Among  the  several  kinds  of 
geometry  design,  the  planar  SOFC  has  higher  energy  conversion 
efficiency  and  a  lower  manufacturing  cost  than  the  tubular  SOFC 
[2,3].  In  order  to  reduce  the  resistance  of  the  ionic  conduction,  an 
extremely  thin  electrolyte  (e.g.  8  pm)  is  usually  employed  in  the 
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positive  electrode,  electrolyte  and  negative  electrode  (PEN)  struc¬ 
ture.  However,  the  thermal  shock  during  the  transient  operations, 
such  as  the  heat-up,  the  start-up,  the  load  change  and  the  cool¬ 
down  processes,  may  damage  the  structure  robustness  of  the 
ultra-thin  electrolyte  and  the  durability  of  the  cell,  since  the  cell 
temperature  varies  in  a  broad  range  from  25  to  800  °C. 

Due  to  the  high  operating  temperature  (700  -  1000  °C),  the 
SOFC  has  to  undergo  the  heat-up  and  start-up  processes  to  warm¬ 
up  the  cell  from  the  room  temperature  (25  °C)  to  the  operable  level 
(around  800  °C).  If  the  cell  temperature  rises  too  fast,  it  will  lead  to 
a  serious  internal  temperature-gradient  increase  and  significant 
differences  of  the  thermal-expansion  rates  of  the  cell  components, 
resulting  in  a  high  thermal  stress  [4]  and  reducing  the  cell  perfor¬ 
mance  and  its  durability.  As  this  is  especially  true  for  the  planar 
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Nomenclature 

Patm 

Atmospheric  pressure  (atm) 

C2 

Inertial  resistance  factor  (m-1) 

P?h 

Partial  pressure  for  the  ith  species  near  the  triple  phase 
boundary  (atm) 

Cp 

Specific  heat  (J  kg-1  K-1) 

Q 

Thermal  energy  (J) 

Dt 

Diffusivity  for  the  ith  species  (m2  s-1) 

Ri 

Chemical  reaction  rate  (mol  m-3  s-1) 

DP 

Particle  diameter  in  porous  zone  (pm) 

i^ohm 

Ohmic  resistance  (Q  cm-2) 

De  ff 

Effective  diffusivity  (m2  s-1) 

R 

Gas  constant  (J  mol-1  K-1) 

E 

Total  energy  (J) 

ofcpl?) 

Fluid  enthalpy  source  term  (J  m-3  s-1) 

Fn 

Theoretical  potential  (V) 

Activation  energy  (kj  mol-1) 

+  V  •  (eplf~v)  =  — rVp  +  V-  (rf)  +  SM  Standard  state 

entropy  (J  mol-1  I< 

Ea 

at 

^anode 

Activation  energy  for  exchange  current  density  of 

-1) 

anode  (kj  mol-1) 

5m 

Momentum-sink  term  in  porous  (N  m  3) 

^cathode 

Activation  energy  for  exchange  current  density  of 

T 

Temperature  (K) 

cathode  (kj  mol-1) 

X 

X  coordinate  (mm) 

Ef 

Total  energy  of  fluid  (J) 

X 

x/L,  dimensionless  coordinate  along  the  cell  length  (-) 

£ocv 

Open  circuit  voltage  (V) 

^cell 

Cell  voltage  (V) 

F 

Faraday  constant  (C  mol-1) 

Y, 

Mass  fraction  for  the  ith  species 

G° 

Gibbs  free  energy  for  the  electrochemical  reaction 

0  mor1) 

Greek  letters 

Hox 

Enthalpy  of  oxidation  (J  mol-1) 

a 

Permeability  (m2) 

hi 

Species  enthalpy  (J  mol-1) 

p 

Transfer  coefficient 

hj 

i 

Standard  state  enthalpy  (J  mol-1) 

T 

The  extent  of  the  anode-recycling 

Current  density  (A  cm-2) 

£ 

Porosity 

io 

Exchange  current  density  (A  cm-2) 

Pact 

Activation  overpotential  (V) 

ji 

Mass  diffusion  flux  for  ith  species  (kg  m-2  s-1) 

Pconc 

Concentration  overpotential  (V) 

k 

Rate  constant  (mol  m-2  bar-1  s-1) 

F 

Dynamic  viscosity  (Pa  s) 

/<o 

Pre-exponential  factor  (mol  m-2  bar-1  s-1) 

p 

Stoichiometric  coefficient  for  j  reactant 

ka 

Pre-exponential  factor  for  exchange  current  density  of 

Stoichiometric  coefficient  for  j  product 

anode  (A  cm-2) 

Density  (kg  m-3) 

kc 

Pre-exponential  factor  for  exchange  current  density  of 

f 

Stress  tensor  (-) 

cathode  (A  cm-2) 

V 

Velocity  (m  s-1) 

ke  ff 

kf 

Effective  thermal  conductivity  (W  m-1  I<-1) 

Fluid  thermal  conductivity  (W  m-1  K-1) 

Abbreviation 

k  s 

Solid  thermal  conductivity  (W  m-1  K-1) 

APU 

Auxiliary  power  unit 

f^WGS 

Equilibrium  constant  for  water  gas  shift  reaction 

DISR 

Direct  internal  steam  reforming 

L 

Cell  length  (mm) 

EC 

Electrochemical 

fn 

Mass  flow  rate  (kg  s-1) 

LSM 

Lanthanum  strontium  manganite 

Mi 

Molecular  weight  (g  mol-1) 

PEN 

Positive-electrolyte-negative 

n 

Number  of  electron  transfer 

TPB 

Triple-phase-boundary 

N 

Number  of  species  (-) 

WGS 

Water  gas  shift  reaction 

P 

Pressure  (atm) 

YSZ 

Yttria-stablized  zirconia 

SOFC,  the  management  of  the  temperature  uniformity  of  the  planar 
SOFC  deserves  in-depth  investigations. 

The  start-up  process  is  usually  initiated  by  feeding  the  fuel  into 
the  fuel  channel  after  the  SOFC  is  heated  up  to  around  600  °C  by 
a  proper  heating  strategy  as  investigated  in  our  previous  work  [5]. 
There  are  two  kinds  of  fuel  reported  to  be  employed  in  the  start-up 
process,  i.e.,  hydrogen  and  hydrocarbon  [6-10].  For  instance,  Seli- 
movic  et  al.  [9]  investigated  the  start-up  time  and  the  maximum 
temperature  gradient  by  feeding  hydrogen  and  methane  as  the  fuel, 
respectively,  for  the  start-up  process  of  a  planar  SOFC.  It  was  found 
that  the  required  time  for  the  case  using  methane  (1  h)  is  2-fold 
longer  than  that  using  hydrogen  (30  min).  For  the  case  using 
hydrogen,  the  maximum  temperature  gradient  occurs  at  the 
beginning  of  the  start-up  process,  with  the  value  about  5  K  cm"1, 
and  drastically  drops  thereafter.  For  the  case  using  methane,  the 
temperature  gradient  keeps  increasing  during  the  whole  start-up 
process,  and  reaches  the  maximum  value,  about  11  I<  cm-1  at  the 
end  of  the  start-up  process.  Petruzzi  et  al.  [8]  showed  that  the 
temperature  of  600  °C  is  the  most  appropriate  state  to  proceed  the 
start-up  process.  In  their  3.5  kW  auxiliary  power  unit  (APU)  system, 


the  start-up  process  takes  about  35  min,  which  is  too  slow  for  their 
APU  system.  Lin  and  Hong  [11]  built  a  250  kW  turbo-SOFC-system 
model  to  investigate  the  start-up  process.  As  the  fuel  inlet  pressure 
is  2  atm,  the  start-up  process  takes  about  1.3  h.  They  also  showed 
that  the  most  important  factor  affecting  the  start-up  time  is  the 
response  of  the  stack  temperature  of  the  SOFC  system.  They  also 
suggested  a  practically  applicable  method  to  accelerate  the  start-up 
process,  i.e.,  increasing  the  inlet  pressure.  As  the  inlet  pressure 
increases  from  2  atm  to  3  atm,  the  start-up  time  is  shortened  to  1  h. 
However,  the  temperature-gradient  variations  are  not  investigated 
in  their  study.  Lin  and  Hong  [12]  investigated  the  cold  start-up 
process  of  a  100  W  SOFC-APU  system  by  their  in-house  devel¬ 
oped  numerical  model.  In  their  system,  the  embedded  burner  and 
the  recycling  mechanism  were  incorporated  in  their  model.  By 
neglecting  the  effect  of  the  temperature  gradient,  it  takes  100  s  for 
the  SOFC-APU  system  to  reach  the  steady  state.  However,  they  also 
mentioned  that  such  a  fast  start-up  pace  will  lead  to  the  problems 
of  the  significant  temperature  gradient,  thermal  stress,  and  fatal 
component  fracture.  Therefore,  a  nonlinear  sliding-observer  was 
developed  in  their  model  to  obtain  the  instantaneous  variation  of 
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the  cell  temperature.  They  also  mentioned  that  such  an  observer  is 
a  useful  tool  for  the  development  of  a  quick  start-up  system.  By 
feeding  back  the  instantaneous  cell  temperature  to  the  control  unit, 
the  high  temperature  gradient  can  be  avoided  in  the  start-up 
process.  Wang  et  al.  [10]  concluded  that  the  dynamic  response 
and  the  structure  integrity  are  significant  issues  for  the  quick  start¬ 
up  process  of  the  SOFC-APU  system.  Cheekatamarla  et  al.  [7] 
experimented  with  the  start-up  process  of  a  SOFC  system  by 
feeding  hydrogen  into  the  cell  at  first.  As  the  cell  temperature 
reached  to  the  normal  operation  level  (~800  °C),  the  fuel  was 
switched  to  propane.  Barzi  et  al.  [6]  used  hydrogen  as  the  fuel  to 
investigate  the  start-up  process  of  a  tubular  SOFC.  With  the 
condition  of  70%  of  the  fuel  utilization  rate  and  800  °C  as  the  inlet 
temperature,  the  start-up  process  takes  about  3  h.  However,  the 
corresponding  temperature  gradient  within  the  cell  was  not  eval¬ 
uated.  In  order  to  reduce  the  temperature  gradient  in  the  cell, 
Nikooyeh  et  al.  [13]  incorporated  the  anode-recycling  mechanism 
in  the  steady-state  analysis  of  the  SOFC  system.  It  was  found  that 
the  temperature  gradient  reduces  about  38%  as  the  anode-recycling 
level  is  75%.  However,  the  effect  of  anode-recycling  mechanism  on 
the  start-up  process  was  not  conducted. 

For  the  above  reviewed  studies  about  the  start-up  process,  it 
was  found  that  the  start-up  pace  is  generally  too  slow,  especially  for 
the  APU  system.  As  the  most  studies  were  focused  on  the  start-up 
time,  the  variation  of  the  temperature  gradient  in  the  SOFC  was  not 
comprehensively  evaluated  to  a  satisfactory  level.  A  complete 
evaluation  for  the  performance  of  the  start-up  process  should  take 
into  account  not  only  the  start-up  time  but  also  the  detailed 
distribution  of  the  temperature  gradient. 

In  the  present  study,  a  finite-volume  model  is  developed  to 
investigate  the  temperature  variation,  the  start-up  time  and  the 
corresponding  temperature  gradient  for  an  anode-supported 
planar  SOFC  during  the  start-up  process.  A  comprehensive 
numerical  simulation  model,  employing  the  finite-volume 
approach,  has  been  developed  for  the  three-dimensional,  multi- 
component,  electrochemically,  and  chemically  reacting  flow  of 
a  single  planar  SOFC  unit  cell.  Three  start-up  procedures  are 
proposed  and  investigated.  The  first  is  to  introduce  the  inlet  fuel  at 
the  operation  temperature  of  800  °C  after  the  heat-up  process  is 
completed.  The  second  is  to  incorporate  the  anode-recycling 
mechanism  into  the  start-up  process.  The  third  is  to  fix  the  differ¬ 
ence  between  the  inlet-fuel  temperature  and  the  cell  minimum 
temperature.  The  effect  of  the  fuel  type  employed  on  the  start-up 
process  is  also  evaluated. 


2.  Model  and  formulation 

2.1.  Start-up  approach 

For  the  start-up  process,  the  steam/carbon  ratio  and  mass  flow 
rate  are  selected  to  keep  the  start-up  time  within  1  h.  To  initiate  the 
start-up  process,  the  fuel  stream  is  fed  into  the  fuel  channel  and  air 
is  fed  into  the  air  channel.  For  the  thermal  condition,  the  inlet 
temperature  is  switched  from  the  state  of  the  end  of  heat-up 
process  to  1073  K  at  the  beginning  of  the  start-up  process  [6]. 
Two  kinds  of  fuel,  i.e.,  hydrogen  and  methane  are  employed  for  the 
start-up  process.  For  the  case  using  hydrogen,  only  the  electro¬ 
chemical  reaction  occurs  within  the  SOFC.  For  the  case  using 
methane,  the  direct  internal  steam  reforming  (DISR),  water  gas 
shift  (WGS)  and  the  electrochemical  reaction  are  considered  in  the 
start-up  process.  The  initial  temperature  distribution  is  adopted 
from  the  optimal  heat-up  results  of  our  previous  study  [5].  The 
start-up  process  finishes  as  the  average  cell  temperature  reaches  to 
99%  of  the  steady-state  value. 


2.2.  Anode  supported  planar  SOFC 

A  finite-volume,  three-dimensional  numerical  model  for  an 
anode-supported  planar  SOFC  cell  unit  is  developed  using  the 
commercial  software  FLUENT  with  the  Gambit  as  the  preprocessor 
[14].  The  grid  system  of  the  multi-zones  for  the  anode  channel,  the 
cathode  channel,  the  interconnector  and  the  PEN  structure,  is 
constructed  by  the  Gambit.  The  conservation  equations  of  mass, 
species,  momentum  and  energy  with  proper  designated  boundary 
conditions  are  solved  by  the  FLUENT. 

The  geometry  of  the  investigated  planar  SOFC  is  shown  in  Fig.  la 
and  the  component  specification  is  shown  in  Fig.  lb.  There  are  20 
unit  cells  on  each  single  cell  plate,  and  10  repeated  single  cell  plate 
for  the  whole  SOFC  stack.  With  the  power  density  of  0.62  W  cm-2, 
the  power  of  the  SOFC  stack  is  about  620  W.  Due  to  the  planar 
geometry,  it  is  commonly  assumed  that  each  unit  cell  has  the  same 
properties,  i.e.,  species  and  temperature  distributions.  It  is  assumed 
that  the  stack  is  well  insulated  [8]  during  the  transient  operations, 
and  the  heat  loss  of  the  unit  cell  at  the  edge  is  negligible.  Due  to  the 
symmetric  boundary  conditions  for  each  cell,  all  the  cells  will  thus 
behave  in  the  same  way  approximately.  A  quantitative  verification 
[13]  also  showed  that  the  temperature  variation  is  within  1  K  for 
over  90%  cells.  Therefore,  one  can  compute  only  one  unit  cell  and 
extend  the  results  to  the  whole  plate  of  the  cell,  as  shown  in  Fig.  la. 
The  geometry  data  of  the  planar  SOFC  is  shown  in  Table  1.  With  the 
symmetric  feature  of  the  planar  SOFC,  only  a  half  of  the  unit  cell  is 
built  in  the  finite-volume  model.  In  the  following  results,  the 
dimension  is  normalized  from  0  to  1.  The  detail  parameters  for  the 
SOFC  system  specification  are  listed  in  Table  1. 


z 


Fig.  1.  Schematic  diagram  of  the  anode-supported  planar  SOFC  for  (a)  stack,  (b)  unit 
cell  and  component  specifications. 
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Table  1 

The  SOFC  system  specifications. 


Anode  thickness 

1000  pm 

Electrolyte  thickness 

10  pm 

Cathode  thickness 

50  pm 

Cell  dimension 

10  cm  x  10  cm 

Outer  channel  dimension 

5  mm  x  2  mm 

Inner  channel  dimension 

4  mm  x  1.5  mm 

Air  to  fuel  ratio 

8 

Steam  to  carbon  ratio 

2 

Fuel  utilization  rate 

70% 

Inlet  temperature 

1073  K 

Operating  Pressure 

1  atm 

2.3.  Conservation  equations 


The  conservation  equations  for  the  start-up  operation  are 
described  in  terms  of  the  porosity.  For  the  non-porous  regions,  the 
porosity  is  equal  to  unity.  The  porosities  for  various  materials  are 
listed  in  Table  2.  The  mass  conservation  equation  for  the  finite 
volume  model  is: 


Q  Op) 

dt 


+  V-  (cpTT)  =  0 


(i) 


The  momentum  conservation  equation  is: 


9  (£p~v) 

dt 


+  V- (cpT^V)  =  -cVp  +  V-  (et)  +  S 


M 


(2) 


where  r  is  stress  tensor,  and  Sm  is  the  momentum  sink  term  for  the 
porous  media.  The  latter  is  evaluated  by: 


Sm  = 


^ ~vi  +  C^pVmzgVi 


(3) 


Cz 


The  inertial  resistance  C2  is  evaluated  by: 

=  3.5  (1  —  e) 

Dn  £3 


(4) 


where  Dp  is  the  average  particle  diameter,  and  e  is  the  porosity. 
Those  parameters  are  listed  in  Table  2. 

For  the  species  conservation  equation,  the  source/sink  term  for 
the  chemical  and  electrochemical  reactions  are  considered  [15]: 

+  V  •  0  pVYi)  =  V  •  0  pDjVV,)  +  jr  MiRi  +  Si  (5) 

1=1 


Table  2 

Physical  properties. 


Property 

Anode 

(Ni-doped 

YSZ) 

Cathode 

(LSM) 

Electrolyte 

(YSZ) 

Interconnector 

(Stainless 

steel) 

Density 

3030  [34] 

3310  [34] 

5163  [34] 

8030  [34] 

(kg  m-3) 

Conductivity 

5.84  [34] 

1.86  [34] 

2.16  [34] 

20.0  [34] 

(W  ITT1  K"1) 

Specific  heat 

595  [34] 

573  [34] 

606  [34] 

502  [34] 

(J  kg-1  K-1) 

Permeability  (-) 

3.4  x  10-2 
[35] 

3.7  x  10  2 

[35]  1.0  x  10  18  [31] 

1.0  x  10~18 

[31] 

Porosity  (-) 

0.42  [35] 

0.36  [35] 

1.0  x  10  5 

[31] 

1.0  x  10~5 

[31] 

Average  particle 

1.40a 

1.40a 

1.0a 

1.0a 

diameter  (pm) 

a  Assumed  value. 


where  Ri  is  the  species  production  by  chemical  reactions,  i.e.,  DISR 
and  WGS,  as  described  above,  and  S;  is  the  sink/source  due  to  H2 
electrochemical  reaction.  The  latter  can  be  expressed  according  to 
the  Faraday’s  law  as: 

Anode  :  %  =  ^  (6) 


(7) 


Cathode :  Sq2 


-i 

nF 


(8) 


where  n  is  the  number  of  electrons  transferred  in  the  electro¬ 
chemical  reaction  with  2  for  the  anode  and  4  for  the  cathode.  The 
energy  conservation  equation  also  takes  the  chemical  and  elec¬ 
trochemical  reactions  into  account,  and  is  given  as  [15]: 


9  (epE) 
dt 


+  V-[  v  ( epE  +  sp)\ 


=  V- 


r/<effVT 


+  (rr-  v ) 


N 


E  £h°Ri + s 


7 


(9) 


where  S5  is  the  fluid  enthalpy  source  term  which  will  be  described 
below.  The  effective  conductivity  in  the  porous  media  is  evaluated 
from  the  fluid  conductivity,  the  solid  conductivity,  and  the  porosity 
as  [16]: 

/<eff  —  £kf  +  (1  —  c)ks  (10) 

In  the  energy  equation,  the  energy  source  terms  due  to  the 
ohmic  heating  corresponding  to  the  electrochemical  reaction  can 
be  expressed  as: 


S 


energy  _ohm 


(ii) 


Another  energy  source  term  in  Eq.  (9)  is  due  to  the  electro¬ 
chemical  reaction  and  can  be  expressed  as  [13]: 


S 


energy -electrochemical 


V  nF 


(12) 


where  A H0x  is  the  enthalpy  of  oxidation.  Since  the  inner  temper¬ 
ature  gradient  is  modest,  the  interior  radiation  is  negligible  [17].  On 
the  outer  surface,  a  properly  designed  insulation  case  is  supposed 
to  keep  the  stack  wall  temperature  below  400  K  [8].  Furthermore, 
the  quantitative  evaluation  showed  that  the  radiation  was  minor  in 
the  SOFC,  specifically,  only  about  0.6%  [18].  Therefore,  the  thermal 
radiation  is  not  considered  in  the  present  model. 


2.4.  Chemical  reactions 

In  the  fuel  side,  the  fast  and  endothermic  DISR  reaction  (Eq.  (13)) 
converts  CF14  to  H2  for  the  electrochemical  reaction  and  CO  for  the 
subsequent  WGS  reaction.  The  latter  reaction  is  to  generate  more 
H2  and  some  CO2  (Eq.  (14)). 

CH4  +  H20~3H2  +  C0  (13) 

C0  +  H20~H2  +  C02  (14) 

The  DISR  reaction  includes  many  intermediate  radicals  and 
complicated  intermediate  reactions,  such  as  absorption,  desorption 
and  decomposition.  In  the  present  study,  the  DISR  reaction  is 
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reduced  to  a  two-equation  mode,  since  it  has  been  proven  to  be 
able  to  provide  adequate  predictions  for  the  simulation  of  the  SOFC 
[11,13,19-21].  The  rate  constants  for  the  above  reactions  are 
described  by  the  Arrhenius  expression: 

k  =  k0e~Ea/RT  (15) 

Although  the  WGS  is  usually  assumed  to  be  equilibrium  [22], 
a  simple  kinetic  model  including  forward  and  backward  reactions  is 
employed  in  the  present  model  in  order  to  account  for  the  dynamic 
change  of  the  species  in  the  start-up  process.  The  related  parame¬ 
ters  for  the  chemical  reactions  of  the  DISR  and  WGS  reactions 
applied  in  the  present  FV  model  are  [23,24]: 


Vir  =  4274  mol  m  2  bar  1  s  1 

(16a) 

Ea,iR  =  82  kj  mol-1 

(16b) 

^0, WGS, forward  =  0.0171  mol  ITT3  Pa~2  S"1 

(17a) 

G, wgs  =  103  kj  mol  1 

(17b) 

The  backward  reaction  rate  of  WGS  is  derived  from  the  equi¬ 
librium  constant  as: 


^0, WGS, backward 


^-0, WGS,  forward 

Kwgs 


The  equilibrium  constant  is  evaluated  by: 


(18) 


^WGS 


exp 


'AS0 

^WGS 

R 


^WGs]  /Patm\  (VJ  vj) 

RT  J\RTJ 


(19) 


where  v'-  is  the  stoichiometric  coefficient  of  the  product  and  v'-  is 
that  of  the  reactant.  The  term  patm  denotes  the  atmospheric  pres¬ 
sure.  The  term  within  the  exponential  function  represents  the 
change  in  Gibbs  free  energy,  and  can  be  computed  as  follows: 


AS0 

^3wgs 

R 


R 


(20) 


A  H° 
ajlWGS 

RT 


N  ,  A  h° 

£('?-■!)  i 

J=  1 


(21) 


where  S?  and  h?  are  the  standard-state  entropy  and  standard-state 
enthalpy,  respectively. 


2.5.  Electrochemical  reactions 

The  hydrogen  generated  via  the  DISR  and  WGS  reactions 
described  above  then  undergoes  the  electrochemical  reaction 
within  the  PEN  structure  as  follows: 

Anode:  H2  +  02~  ->H20  +  2e~  (22) 


The  electrochemical  reaction  is  assumed  to  take  place  near 
the  (TPB)  triple-phase-boundary  within  the  anode  layer  for  the 
present  anode-supported  design.  Basically,  CO  can  also  carry  out 
the  electrochemical  reaction  in  a  SOFC,  but  the  reaction  rate  is 
much  lower  than  that  of  H2  oxidation  [22,25].  Therefore,  only  the 
electrochemical  reaction  of  H2  is  considered  in  the  present 
model. 


2.6.  Electrochemical  model 


The  cell  voltage  is  related  to  the  open  circuit  voltage  and  three 
polarization  loss  terms.  The  open  circuit  voltage  is  calculated  by  the 
Nernst  equation: 


£qcv  = 


c  RTi 

E0  +  -=ln 
nF 


tpb  tpb1/2' 
*  H2  *  O2 
„tpb 

Ph2o 


(25) 


The  open  circuit  voltage  is  affected  by  the  partial  pressure  of 
related  species  in  the  triple  phase  boundary.  The  theoretical 
potential  Eq  is  given  as 


-AG° 

nF 


(26) 


where  the  AG°  is  the  Gibb’s  free  energy  for  the  electrochemical 
reaction  of  hydrogen.  The  cell  voltage  is  evaluated  by: 


^cell  —  Rr 


ocv  lRohm  Vacfa  Vacfc  Vconc 


(27) 


where  ^act,a  and  ^actjC  stand  for  the  activation  loss  for  the  anode  and 
the  cathode,  respectively.  It  can  be  described  by  the  Butler- Volmer 
equation: 


*cell  —  *0,  electrode 


exp 


7 SnFrj 


act, electrode 


RT 


-  exp 


(1  @)nFriac  t, electrode  \ 

Or  )\ 


(28) 


where  the  electrode  includes  the  anode  and  the  cathode.  There 
were  two  approaches  reported  to  estimate  the  activation  polar¬ 
ization  [26].  For  the  first  approach,  the  expressions  for  activation 
polarizations  of  the  two  electrodes  are  different  [27].  For  the 
second  approach,  the  transfer  coefficient  (/?  in  Eq.  (28))  is 
assumed  as  0.5  [13,22,28-32].  Therefore,  the  Butler-Volmer 
equation  (Eq.  (28))  can  be  rewritten  into  an  explicit  form,  and 
the  same  expression  can  be  employed  to  both  electrodes  with 
the  different  exchange  current  densities.  In  the  present  study, 
the  second  approach  is  employed  in  the  electrochemical  model 
and  the  exchange  current  density  in  Eq.  (28)  is  determined 
by  [29]: 


,  /PTH?UPS?o\c„D/-£act,anode\  .  . 

*0, anode  —  ^anode  I  p  J  [  p  I  CXp)  (  ^  J  (2)3) 

(Tpb\  0-25 

P02  \  f—^act, cathode \ 

p  j  (  RT  J  (30) 


1  9 

Cathode:  -02  +  2e  ^O2 

2 

Net:  H2+1o2^H20 


where  the  pre-exponential  factor  (K)  is  5.5  x  108  A  m~2  for  the 

(23 )  anode  and  7  x  10s  A  m-2  for  the  cathode  [29].  The  activation  energy 
(Eact)  is  100  x  103  J  mor1  for  the  anode  and  120  x  103  J  mol-1  for 
the  cathode  [29].  The  concentration  loss  in  Eq.  (27)  is  due  to  mass 

(24)  diffusion  for  the  reactants  and  products  within  the  electrode,  and 
can  be  described  by  [29]: 
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RT]n  ( PHj>go\ 

nF  [pgPKo) 


(31) 


Since  the  cell  is  anode-supported,  the  concentration  loss  on  the 
cathode  side  is  negligible  [22].  The  parameters  for  the  activation 
loss  are  adopted  from  the  study  of  Costamagna  et  al.  [29,33].  The 
physical  properties  of  the  materials  used  in  the  present  model  are 
given  in  Table  2. 

For  the  numerical  method  of  the  FV  model,  the  scalars  will  be 
stored  at  the  cell  center,  while  face  values  are  interpolated  from  the 
cell-center  value  using  the  upwind  scheme.  The  upwind  scheme 
interpolation  is  employed  for  the  density,  momentum,  species  and 
energy  equations.  The  computational  grid  was  tested  to  confirm  the 
grid  independency.  The  government  equations  are  solved  by  the 
pressure-based  solver  [14].  Therefore,  the  momentum  and  conti¬ 
nuity  equations  are  solved  simultaneously  as  step  one.  After 
updating  the  mass  flux,  the  species  and  energy  equations  are  solved 
as  step  two.  If  the  convergence  criterion  is  met,  the  iteration  will  be 
terminated.  Otherwise  the  properties  are  updated  and  back  to  step 
one  to  initiate  a  new  iteration  loop  until  the  convergence  criterion 
is  met. 

2.7.  Anode -re cycling  mechanism 

The  anode-recycling  mechanism  is  employed  in  the  start-up 
process  by  mixing  the  anode  exhaust  gas  with  fuel  before  the 
SOFC  module  [1  ].  It  is  expected  to  reduce  the  required  start-up  time 
and  the  corresponding  temperature  gradient.  The  extent  of  the 
anode-recycling  is  defined  as: 

Yi,y%  =  ^i',out  x  y%  +  YI;inlet  x  0  —  7%)  (32) 


effect  on  the  start-up  process.  Finally,  the  mechanism  of  the  fixed- 
temperature-difference  is  employed  to  explore  whether  it  is  able  to 
reduce  the  start-up  time  and  the  effective  maximum  absolute 
temperature-gradient  (the  largest  one  during  the  whole  process). 
The  start-up  time  and  the  effective  maximum  absolute 
temperature-gradient  are  taken  as  the  criteria  to  evaluate  the 
performance  of  the  start-up  mechanism. 


3.  Results  and  discussion 

3.1.  The  validation  of  the  finite-volume  model 

Numerical  validation  is  made  for  the  FV  model.  The  calculated 
current-voltage  curve  is  first  compared  with  the  experiment  data 
[29],  as  shown  in  Fig.  2.  The  trends  of  the  predicted  and  experiment 
results  match  with  each  other.  The  difference  is  about  3%  within  the 
range  from  0.5  to  1.0  V,  validating  the  electrochemical  model.  Next, 
the  temperature  profile,  as  shown  in  Fig.  3,  is  compared  with  that  of 
Nikooyeh  et  al.  [  13  ]  for  the  case  of  V  =  0.7  V  to  validate  the  chemical 
reaction  for  methane  feedstock.  The  inlet  conditions,  geometric 
parameters,  and  operating  conditions  are  adopted  from  the 
compared  reference  [13].  Since  the  material  properties  were  not 
reported  in  the  paper  of  Nikooyeh  et  al.  [13],  typical  material 
properties  for  the  anode  (Ni-YSZ),  the  electrolyte  (YSZ)  and  the 
cathode  (LSM)  are  assumed  in  the  present  study  [31,34,35].  The 
interconnector  uses  stainless  steel  [34].  As  shown  in  Fig.  3,  the 
temperature  profiles  predicted  in  the  present  study  are  in  good 
agreement  with  those  by  Nikooyeh  et  al.  [13].  The  overall  temper¬ 
ature  difference  is  less  than  1%.  Furthermore,  the  maximum 
temperature  gradient  (17  K  cm-1),  shown  in  Fig.  3,  also  matches 
closely  with  the  compared  data  (16  K  cm-1)  [13],  validating  the 
present  FV  model. 


Under  the  anode-recycling  condition,  a  specified  amount,  i.e.,  y 
%  in  Eq.  (32),  of  the  inlet  species  is  replaced  by  the  outlet  gas.  The 
anode-recycling  mechanism  is  first  employed  to  the  steady-state 
simulation  and  compared  with  the  reference  result  as  a  valida¬ 
tion.  Next,  it  will  be  incorporated  into  the  start-up  process  for  the 
case  using  methane. 

2.8.  Fixed-temperature-difference  mechanism 

For  the  fixed-temperature-difference  mechanism,  the  temper¬ 
ature  difference  to  be  fixed  is  defined  as  the  difference  between  the 
inlet  temperature  and  the  minimum  temperature  of  the  PEN 
structure  as: 

^inlet  =  ^celLmin  ~F  AT  (33) 

This  mechanism  is  employed  to  the  case  utilizing  the  hydrogen 
and  methane  to  examine  its  performance  on  the  start-up  process. 

2.9.  Simulation  approaches  and  procedures 

In  the  present  study,  the  results  of  the  thermal  behavior,  the 
start-up  time  and  the  corresponding  temperature  gradient  are 
investigated  for  an  anode-supported  planar  SOFC  during  the  start¬ 
up  process.  The  validation  of  the  proposed  model  is  presented  by 
comparing  the  polarization  curve  of  the  cell  and  the  cell  temper¬ 
ature  profile  with  those  of  relevant  studies.  The  hydrogen  and 
methane  are  employed  in  the  start-up  process  of  the  SOFC.  First,  the 
effect  of  different  fuels  on  the  start-up  performance  is  investigated. 
Next,  the  anode-recycling  mechanism  is  incorporated  to  verity  its 


3.2.  Effect  of  different  fuel  on  the  start-up  process 

The  effect  of  different  fuels,  i.e.,  hydrogen  and  methane,  on  the 
start-up  process  is  investigated  in  this  section.  For  the  case  using 
hydrogen,  the  variation  of  the  average  cell-temperature  during  the 
start-up  process  is  shown  in  Fig.  4a,  while  the  variation  of  the 
effective  maximum  absolute  temperature-gradient  for  the  first  10  s 
is  shown  in  Fig.  4b.  At  the  beginning  of  the  start-up  process,  the 
average  cell-temperature  increases  relatively  fast  (Fig.  4a),  leading 
to  the  largest  effective  maximum  absolute  temperature-gradient 
(Fig.  4b).  When  the  increasing  cell-temperature  is  close  to  the 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1  1.1  1.2 

Current  Density  (A  cm  2) 

Fig.  2.  Comparison  of  the  predicted  current-voltage  curve  with  the  experimental 
data  [29]. 
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Fig.  3.  Temperature  profiles  under  various  operating  conditions  and  the  corresponding 
temperature-gradient  distributions  in  comparison  with  Nikooyeh  et  al.  [13]. 


inlet  temperature  (1073  K),  the  start-up  pace  slows  down  and 
a  smaller  effective  maximum  absolute  temperature-gradient  is 
exhibited  (Fig.  4b),  since  the  difference  between  the  cell  tempera¬ 
ture  and  the  inlet  temperature  reduces.  The  whole  start-up  process 
takes  about  14  min.  The  largest  effective  maximum  absolute 
temperature-gradient  is  about  113.6  K  cm-1,  which  occurs  at  about 
1.6  s  from  the  beginning  of  the  start-up  process.  The  temporal 
variation  of  the  cell-temperature  distribution  during  the  start-up 
process  is  shown  in  Fig.  5a,  while  the  corresponding  variation  of 
the  temperature-gradient  distribution  is  shown  in  Fig.  5b.  The 
initial  temperature  distribution  (at  0  s)  is  adopted  from  the  heat-up 
results  of  the  counter-flow  configuration  in  our  previous  study  [5], 
where  the  largest  maximum  absolute  temperature-gradient  is 
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Fig.  4.  (a)  Temperature  variation  and  (b)  the  corresponding  variation  of  the  effective 
maximum  absolute  temperature-gradient  within  the  first  10  s  for  the  case  utilizing  the 
hydrogen  in  the  start-up  process. 


Fig.  5.  (a)  Variation  of  the  temperature  distribution  and  (b)  the  corresponding  varia¬ 
tion  of  the  distribution  of  the  temperature  gradient  for  the  case  utilizing  the  hydrogen 
in  the  start-up  process. 


about  16  I<  cm'1  (Fig.  5b).  As  the  high-temperature  hydrogen  is 
introduced  from  the  fuel  inlet,  the  temperature  near  the  inlet 
increases  drastically.  Since  the  downstream  temperature  has  not 
caught  up  with  the  pace  of  the  inlet  region,  the  largest  effective 
maximum  absolute  temperature-gradient  (113.6  I<  cm'1)  is 
exhibited  near  the  inlet  region  at  the  beginning  of  the  start-up 
process  (Fig.  5b).  After  that,  as  the  downstream  temperature 
increases,  the  absolute  temperature-gradient  near  the  inlet  region 
gradually  decreases.  At  about  200  s  from  the  beginning  of  the  start¬ 
up  process,  the  absolute  temperature-gradient  (15  I<  cm-1)  reduces 
to  the  level  of  the  initial  condition  of  the  start-up  process.  At  about 
400  s  from  the  beginning  of  the  start-up  process,  the  downstream 
temperature  gradually  becomes  larger  than  that  near  the  inlet 
region  (Fig.  5a).  At  800  s  from  the  beginning  of  the  start-up  process, 
the  temperature  gradient  turns  to  positive  near  the  inlet  region 
(Fig.  5b)  as  the  downstream  temperature  becomes  larger  than  that 
of  the  upstream  region  (Fig.  5a).  It  is  noticeable  that  the  feeding  fuel 
heats  the  cell  at  the  beginning  of  the  start-up  process,  while  it  cools 
the  cell  at  the  later  stage  of  the  start-up  process.  This  role- 
transformation  was  also  observed  in  the  study  of  Selimovic  et  al. 
[9]  for  the  start-up  process  of  a  SOFC. 

For  the  case  using  methane  to  start-up  the  SOFC,  the  variation  of 
the  average  cell  temperature  during  the  start-up  process  is  shown 
in  Fig.  6a,  while  the  variation  of  the  effective  maximum  absolute 
temperature-gradient  for  the  first  50  s  is  shown  in  Fig.  6b.  The  trend 
of  the  variation  of  the  average  cell-temperature  is  similar  to  that 
using  hydrogen  with  a  much  slower  pace.  The  whole  process  takes 
about  45  min,  and  the  effective  maximum  absolute  temperature- 
gradient  is  about  53  K  cm-1,  occurring  at  about  9.5  s  from  the 
beginning  of  the  start-up  process.  It  is  noticeable  that  the  average 
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Fig.  6.  Variation  of  the  (a)  temperature  and  (b)  the  corresponding  temperature 
gradient  for  the  case  utilizing  the  methane  in  the  start-up  process. 
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Fig.  7.  (a)  Variation  of  the  temperature  distribution  and  (b)  the  corresponding  varia¬ 
tion  of  the  distribution  of  the  temperature  gradient  for  the  case  utilizing  the  methane 
in  the  start-up  process. 


cell  temperature  drops  slightly  at  the  beginning  of  the  start-up 
process.  This  is  because  the  hydrogen-required,  exothermic  elec¬ 
trochemical  reaction  takes  place  only  after  the  endothermic 
internal  reforming  reaction  has  proceeded,  leading  to  a  small  dip  of 
the  cell  temperature  at  the  beginning  of  the  start-up  process.  The 
temporal  variation  of  the  cell-temperature  distribution  is  shown  in 
Fig.  7a,  while  the  corresponding  variation  of  the  temperature- 
gradient  distribution  is  shown  in  Fig.  7b.  The  temporal  variation 
of  the  temperature  distribution  is  similar  to  that  of  the  case  using 
hydrogen.  In  the  early  stage,  the  temperature  near  the  inlet  region 
increases,  resulting  in  the  largest  effective  maximum  absolute 
temperature-gradient.  However,  the  extent  is  much  smaller  than 
that  of  the  case  utilizing  hydrogen.  The  temperature  rising  pace  is 
slower  for  the  case  using  methane,  since  the  hydrogen  for  the 
exothermic  electrochemical  reaction  is  not  directly  provided,  but 
generated  via  the  internal  reforming  and  the  water-gas-shift  reac¬ 
tion.  Besides,  the  endothermic  internal  reforming  reaction 
moderates  some  extent  of  the  thermal  energy  provided  by  the  gas 
fed  at  the  inlet.  The  effective  maximum  absolute  temperature- 
gradient  for  the  case  utilizing  methane  is  reduced  by  54%, 
compared  to  that  utilizing  hydrogen. 

A  comprehensive  comparison  for  the  effect  of  different  fuels  on 
the  start-up  performance  is  shown  in  Fig.  8.  The  time  required  for 
the  start-up  process  for  the  case  utilizing  methane  is  about  45  min 
and  3.2  times  longer  than  that  using  hydrogen.  This  is  because  the 
hydrogen  for  the  electrochemical  reaction  has  to  be  generated  from 
the  steam-reforming  reaction  for  the  case  utilizing  methane.  For 
the  APU  system,  however,  a  start-up  process  taking  more  than 
40  min  is  too  slow  [8].  The  effective  maximum  absolute 
temperature-gradient  for  the  case  using  hydrogen  is  2.2  times 


larger  than  that  utilizing  methane.  This  is  because  the  endothermic 
internal  reforming  reaction  of  methane  consumes  some  thermal 
energy  from  the  feeding  gas,  narrowing  down  the  temperature 
difference  between  the  inlet  and  the  interior  in  the  early  stage.  At 
the  end  of  the  start-up  process,  the  temperature  gradients  are 
comparable  for  both  fuels.  Therefore  the  effect  of  different  fuels  on 
the  temperature  uniformity  of  the  SOFC  system  is  more  significant 
at  the  beginning  of  the  start-up  process.  The  commonly  referred 
statement  that  the  endothermic  internal  reforming  reaction 
generates  a  large  temperature-gradient  [13,22,24,30,36]  may  be 
true  in  the  steady-state  operation,  while  it  is  a  positive  effect  on  the 


Time  (min)  TGmax  (K  cnr1)  TG  at  end 

(K  cm'1) 

Fig.  8.  Comprehensive  comparison  of  the  results  by  utilizing  different  fuel  in  the  start¬ 
up  process  ( TGmax  is  the  effective  maximum  absolute  temperature-gradient,  and  TG  at 
end  is  the  temperature-gradient  at  the  end  of  the  start-up  process). 
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accommodation  of  the  temperature  uniformity  during  the  start-up 
process.  Therefore,  the  internal  reforming  reaction  is  favorable  in 
the  start-up  process,  since  it  moderates  the  large  temperature 
gradient. 

3.3.  Effect  of  the  anode-cycling  on  the  start-up  process 

The  temperature  distributions  at  the  steady  state  under  the 
anode-recycling  conditions  are  shown  in  Fig.  9a,  while  the 
temperature-gradient  distributions  are  shown  in  Fig.  9b.  The 
temperature  near  the  inlet  region  increases  with  the  increasing 
extent  of  the  anode-recycling  (Fig.  9a).  This  is  because  the  mass 
fraction  of  the  methane  reduces,  leading  to  a  less  extent  of  the 
endothermic  internal  reforming  reaction.  For  the  downstream 
region,  the  temperature  increases  with  the  increasing  extent  of  the 
anode-recycling,  since  some  hydrogen  is  introduced  via  the  recy¬ 
cling,  leading  to  a  larger  current  density  (0.57  A  citT2  without 
anode-recycling  and  0.71  A  cm-2  for  70%  anode-recycling).  For  the 
temperature  gradient,  the  maximum-temperature-gradient 
reduces  with  the  increasing  extent  of  the  anode-recycling,  and  its 
position  shifts  toward  the  inlet  region  (Fig.  9b).  As  the  extent  of  the 
anode-recycling  increases  to  70%,  the  negative  temperature 
gradient  near  the  inlet  region  vanishes.  This  is  because  the  endo¬ 
thermic  internal  reforming  reaction  reduces  with  the  reduction  of 
the  mass  fraction  of  methane.  As  70%  of  the  anode-recycling 
mechanism  is  employed,  the  maximum-temperature-gradient 
reduces  from  17  I<  cm-1  to  12  I<  cm-1  (decrease  5  I<  cm-1),  while 
it  reduces  from  16  K  cm-1  to  10  K  cm-1  (decrease  6  K  cm'1)  in  the 
study  of  Nikooyeh  et  al.  [13].  Therefore,  the  effect  of  the  anode¬ 
recycling  mechanism  at  the  steady  state  proposed  by  the  present 
study  is  consistent  with  that  in  the  literature  [13]. 


0  0.2  0.4  0.6  0.8  1 

x  -  Dimensionless  Distance  (x/L) 


The  anode-recycling  mechanism  is  then  incorporated  into  the 
start-up  process,  and  the  results  about  the  start-up  time  and  the 
effective  maximum  absolute  temperature-gradient  are  shown  in 
Fig.  10.  The  effect  of  the  anode-recycling  on  the  start-up  time  is 
significant.  Compared  to  the  base  case  (without  anode-recycling), 
the  start-up  time  is  shortened  by  29.29%  for  the  30%  recycling, 
40.75%  for  the  50%  recycling,  and  48.58%  for  the  70%  recycling.  This 
is  because  less  methane  and  more  hydrogen  are  introduced  under 
the  anode-recycling  condition.  Therefore,  with  the  anode-recycling 
mechanism,  less  methane  undergoes  the  complicated  reactions  of 
IR,  WGS  and  EC  (Electrochemical),  while  more  hydrogen  is  directly 
converted  to  the  electric  power,  leading  to  a  faster  start-up. 
However,  the  effect  of  the  anode-recycling  mechanism  on  the 
effective  maximum  absolute  temperature-gradient  is  insignificant 
(<0.7%)  in  the  start-up  process.  Therefore,  the  incorporation  of  the 
anode-recycling  mechanism  reduces  the  temperature  gradient  at 
the  steady  state,  while  it  shortens  the  time  required  for  the  start-up 
process. 


3.4.  The  fixed-temperature-difference  start-up  technique 

As  shown  in  the  previous  section,  the  high  temperature- 
gradient,  which  occurs  at  the  beginning  of  the  start-up  process,  is 
due  to  the  sudden  increase  of  the  inlet  temperature  to  the  opera¬ 
tion  value.  An  alternative  start-up  procedure  is  introduced  in  the 
present  section,  i.e.,  the  fixed-temperature-difference  mode.  This 
method  is  expected  to  reduce  the  large  temperature  gradient  at  the 
beginning  of  the  start-up  process,  and  speed  up  the  pace  at  the  later 
stage.  For  the  case  using  hydrogen,  the  temporal  variation  of  the 
average  cell-temperature  and  the  effective  maximum  absolute 
temperature-gradient  under  different  fixed  temperature  differ¬ 
ences  are  shown  in  Fig.  11.  It  is  found  that  for  the  case  of  a  fixed 
100  K  temperature-difference,  the  temperature  rising  rate  is  slower 
than  that  of  the  original  case,  while  the  effective  maximum  abso¬ 
lute  temperature-gradient  is  49.3  K  cm-1  which  is  also  smaller  than 
that  of  the  original  case.  For  the  case  of  a  fixed  175  I<  temperature- 
difference,  the  effective  maximum  absolute  temperature-gradient 
is  comparable  with  that  of  the  original  case.  The  effective 
maximum  absolute  temperature-gradient  for  the  case  of  a  fixed 
250  K  temperature-difference  is,  however,  larger  than  that  of  the 
original  case.  A  comprehensive  comparison  of  the  effective 
maximum  absolute  temperature-gradient  and  the  start-up  time  is 
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Fig.  9.  (a)  Temperature  distributions  and  (b)  temperature  gradient  distributions  under 
the  anode-recycling  condition. 


Fig.  10.  Effect  of  the  anode-recycling  mechanism  on  the  start-up  time  and  the  effective 
maximum  absolute  temperature-gradient. 
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Fig.  11.  Temperature  variation  and  the  corresponding  variation  of  the  effective 
maximum  absolute  temperature-gradient  within  the  first  10  s  for  the  case  utilizing  the 
hydrogen  and  the  fixed  temperature  difference  mode  in  the  start-up  process. 
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Fig.  13.  Temperature  variation  and  the  corresponding  variation  of  the  effective 
maximum  absolute  temperature-gradient  within  the  first  50  s  for  the  case  utilizing  the 
methane  and  the  fixed  temperature  difference  mode  in  the  start-up  process. 


shown  in  Fig.  12.  Comparing  the  case  of  a  fixed  100  I<  temperature- 
difference  with  the  original  one,  it  is  noticeable  that  the  time 
required  for  the  start-up  process  is  comparable,  while  the  effective 
maximum  absolute  temperature-gradient  is  reduced  by  56%.  For 
the  case  of  a  fixed  175  K  temperature-difference,  the  effective 
maximum  absolute  temperature-gradient  is  comparable  with  that 
of  the  original  case,  while  the  start-up  time  is  shortened  by  42%.  For 
the  case  of  a  fixed  250  I<  temperature-difference,  the  effective 
maximum  absolute  temperature-gradient  is  41%  larger  than  that  of 
original  case,  while  the  time  required  for  the  start-up  is  58%  shorter 
than  that  of  the  original  case.  Therefore,  a  properly  selected 
temperature-difference  of  the  inlet  fuel  and  the  cell  may  lead  to 
a  faster  start-up  pace  and  a  smaller  temperature  gradient. 

The  results  for  the  case  using  methane  as  the  fuel  under  the 
mechanism  of  a  fixed  temperature-difference  are  shown  in  Fig.  13. 
It  is  found  that  for  the  case  of  a  fixed  100  K  temperature-difference, 
the  temperature  rising  rate  within  the  first  50  s  is  much  slower  than 
that  of  the  original  case,  while  the  effective  maximum  absolute 
temperature-gradient  is  about  25  K  cm-1  smaller  than  that  of  the 
original  case.  For  the  case  of  a  fixed  175  I<  temperature-difference, 
the  effective  maximum  absolute  temperature-gradient  is  compa¬ 
rable  with  that  of  the  original  case.  The  effective  maximum  abso¬ 
lute  temperature-gradient  for  the  case  of  a  fixed  250  K 
temperature-difference  is,  however,  larger  than  that  of  the  original 
case.  A  comprehensive  comparison  of  the  effective  maximum 
absolute  temperature-gradient  and  the  start-up  time  is  shown  in 
Fig.  14.  For  the  case  of  a  fixed  100  I<  temperature-difference,  the 


time  required  for  the  start-up  process  is  about  37%  and  the  effective 
maximum  absolute  temperature-gradient  about  52%  less  than 
those  of  the  original  case.  For  the  case  of  a  fixed  175  K  temperature- 
difference,  the  effective  maximum  absolute  temperature-gradient 
is  comparable,  while  the  start-up  time  is  reduced  by  57%.  For  the 
case  of  a  fixed  250  K  temperature-difference,  the  effective 
maximum  absolute  temperature-gradient  is  39%  larger  than  that  of 
the  original  case,  while  the  start-up  time  is  reduced  by  67%. 
Therefore,  for  the  case  using  methane,  it  is  interesting  to  note  that 
the  start-up  procedure  of  a  fixed  100  K  temperature-difference  may 
reduce  both  the  start-up  time  and  the  effective  maximum  absolute 
temperature-gradient.  The  reduction  of  the  effective  maximum 
absolute  temperature-gradient  is  due  to  a  relatively  smaller  inlet 
temperature,  which  is  the  main  factor  causing  the  increase  of  the 
absolute  temperature-gradient  at  the  beginning  of  the  start-up 
process.  At  the  later  stage  of  the  start-up  process,  the  fixed- 
temperature-difference  mechanism  accelerates  the  start-up 
process  and  shortens  the  required  time.  In  contrary,  for  the  orig¬ 
inal  start-up  procedure,  the  temperature  difference  reduces 
substantially  as  the  cell  temperature  approaches  the  inlet  temper¬ 
ature,  slowing  down  the  start-up  process.  Therefore,  for  the  start¬ 
up  procedure  using  a  fixed-temperature-difference  mode,  an 
optimal  value  for  the  temperature  difference  could  be  selected  to 
speed  up  the  start-up  process,  while  keeping  the  effective 
maximum  absolute  temperature-gradient  under  the  allowable 
threshold. 
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Fig.  12.  Comprehensive  comparison  of  the  start-up  time  and  the  effective  maximum 
absolute  temperature-gradient  for  the  case  using  hydrogen  for  the  effect  of  the  fixed 
temperature  difference  mechanism. 


Fig.  14.  Comprehensive  comparison  of  the  start-up  time  and  the  effective  maximum 
absolute  temperature-gradient  for  the  case  using  methane  for  the  effect  of  the  fixed 
temperature  difference  mechanism. 
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4.  Conclusions 

For  the  start-up  process  of  a  SOFC  system,  both  methane  and 
hydrogen  utilized  as  the  fuel,  respectively,  have  been  investigated  for 
their  effects  on  the  start-up  process.  The  anode-recycling  mechanism 
is  tested  at  the  steady  state  as  a  validation,  and  then  incorporated  into 
the  start-up  process  to  examine  its  performance.  Finally,  the  perfor¬ 
mance  of  the  fixed-temperature-difference  mechanism  is  evaluated 
for  the  start-up  process.  Some  conclusions  are  drawn  from  the  present 
study  as  follows: 

(1)  For  the  start-up  process  using  hydrogen  as  the  fuel,  the  gas 
introduced  at  the  inlet  warms  up  the  cell  at  the  beginning  of 
the  start-up  process,  and  turns  to  cool  the  cell  at  the  later  stage 
of  the  start-up  process  when  the  temperature  of  the  cell  is 
higher  than  that  of  the  inlet  gas.  For  the  start-up  process  using 
methane  as  the  fuel,  there  is  a  slightly  temperature  dip  at  the 
beginning  of  the  start-up  process  due  to  the  proceeding  of  the 
internal  reforming  reaction.  The  effective  maximum  absolute 
temperature-gradient  is  generated  at  the  beginning  of  the 
start-up  process  for  both  cases.  For  the  present  investigated 
SOFC  configuration,  the  start-up  time  for  the  case  utilizing 
methane  is  3.2  times  longer  than  that  using  hydrogen.  The 
effective  maximum  absolute  temperature-gradient  for  the  case 
using  hydrogen  is  2.2-fold  larger  than  that  utilizing  methane. 
The  endothermic  internal  reforming  reaction  of  methane  has 
a  positive  effect  on  the  accommodation  of  the  temperature 
uniformity  during  the  start-up  process. 

(2 )  The  effect  of  the  anode-recycling  on  the  start-up  time  is  significant. 
Comparing  to  the  base  case,  i.e.,  without  anode-recycling,  the 
start-up  time  is  shortened  by  29.29%  for  the  30%  recycling,  40.75% 
for  the  50%  recycling,  and  48.58%  for  the  70%  anode-recycling  for 
the  present  investigated  SOFC  configuration.  However,  the  effect 
of  the  anode-recycling  mechanism  on  the  effective  maximum 
absolute  temperature-gradient  is  not  significant. 

(3)  For  the  fixed-temperature-difference  mechanism  in  the  start¬ 
up  process,  a  properly  selected  temperature-difference  may 
lead  to  a  shorter  start-up  time  and  a  smaller  effective 
maximum  absolute  temperature-gradient.  The  reduction  of  the 
effective  maximum  absolute  temperature-gradient  is  due  to 
a  smaller  inlet  temperature  in  the  early  stage  of  the  start-up 
process.  In  the  later  stage  of  start-up  process,  the  fixed 
temperature-difference  accelerates  the  process  and  shortens 
the  required  time.  Therefore,  the  start-up  process  can  be 
accelerated  by  choosing  an  optimal  value  for  the  temperature 
difference,  while  keeping  the  effective  maximum  absolute 
temperature-gradient  under  an  allowable  threshold. 
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